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theory 

D. V. SenthUkumar 1 ' 2 J. Kurths 2 ' 3 Q and M. Lakshmanan^] 

1 Centre for Dynamics of Complex Systems, 
14469 Potsdam Germany 

2 Potsdam Institute for Climate Impact Research, 
H413 Potsdam Germany 

3 Institute for Physics, Humboldt University, 
12489 Berlin, Germany 
and 

4 Centre for Nonlinear Dynamics, 
Department of Physics, Bharathidasan University, 
Tiruchirapalli - 620 024, India 

(Dated: June 1, 2009) 

Stability of synchronization in unidirectionally coupled time-delay systems is studied using the 
Krasovskii-Lyapunov theory. We have shown that the same general stability condition is valid for 
different cases, even for the general situation (but with a constraint) where all the coefficients of the 
error equation corresponding to the synchronization manifold are time-dependent. These analytical 
results are also confirmed by numerical simulation of paradigmatic examples. 
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Synchronization of coupled chaotic dynamical systems 
is an active area of research in different branches of sci- 
ence and technology [l], 0] • Different kinds of chaos syn- 
chronization have been reported both theoretically and 
experimentally since its discovery in coupled chaotic os- 
cillators |3(. Recent studies on synchronization have 
been focused on coupled time-delay systems with or 
without time-delay coupling because of its intrinsic na- 
ture of generating high dimensional chaotic signals j4| 
and the ease of experimental realization of these sys- 
tems These systems have potential applications in 
secure communication, cryptography, controlling, long 
term prediction, optimization of nonlinear system per- 
formance, modelling brain activity, pattern recognition, 

etc. (cf. H H & 0, & H Specifically ' ^P^- 

ing synchronization, communication with chaotic lasers 

was demonstrated in Q and digital information at gi- 
gabyte rates was transferred successfully [|[ exploit- 
ing time-delayed feedback to generate high-dimensional, 
high-capacity waveforms at high bandwidths. With this 
kind of applications of high-dimensional chaotic signals 
of time-delay systems, it becomes extremely important 
to establish conditions under which synchronized states 
are asymptotically stable in coupled time-delay systems. 
For this purpose, the Krasvoskii-Lyapunov theory has 
become an extremely useful tool 0, |l(J, EH, (3 E3] ■ 

On the other hand, recently, it has been pointed out 



that the Krasovskii-Lyapunov theory is not suitable for a 
more general case where the error equation corresponding 
to the synchronization manifold is time-dependent; espe- 
cially when all its coefficients are time-dependent [111 ]. 
In this contribution we will show that the Krasovskii- 
Lyapunov theory is not restricted to rather special cases 
discussed in the literature so far, but that it can be ex- 
ploited as a powerful tool in identifying the synchroniza- 
tion thresholds and the stability of synchronization in 
general coupled time-delay systems. In particular, we 
will show that the same general stability condition re- 
sulting from the Krasovskii-Lyapunov theory is indeed 
valid for the general case, where the coefficients are time 
dependent. This general situation is very important in 
many ap plic ations such as synchronization via dynamical 
relaying [12j |. for instance, synchronization of neural ac- 
tivity in brain with complex functional architecture has 
been shown to underlie cognitive acts [HI , in dynamically 
evolving networks [l4[ such as ad- hoc networks [l5l |. etc. 
In particular, we will discuss all the four possible cases 
that arise due to the nature of the coefficients in the error 
equation corresponding to the synchronization manifold 
and show that the same stability condition deduced from 
the Krasovskii-Lyapunov functional approach is valid for 
all the cases, subject to certain conditions. We will also 
confirm these analytical results by numerical analysis us- 
ing paradigmatic examples. 

Consider the following linearly coupled scalar time- 
delay system, 
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where a and b are positive constants, r > is the delay- 
time, K(t) is the coupling function between the drive and 
the response systems and f(x) is some nonlinear function. 
Now we can deduce the stability condition for complete 
synchronization of the general unidircctionally coupled 
time-delay systems (p}. The time evolution of the dif- 
ference system with the state variable A = x(t) — y(t) 
(the error equation corresponding to the complete syn- 
chronization manifold of the coupled time-delay system 
(Q])) for small values of A can be written as 

A = — (a + K(t))A + bf\y(t - t))A t , A t = A(f - r). 

(2) 

It is to be noted that there arises four cases depending 
on the nature of the coefficient of the A and A r terms of 
the above error equation as follows: 

1. Both coefficients of the A and A T terms are time- 
independent. 

2. The coefficient of the A term is time-independent 
and that of the A T term is time-dependent. 

3. The coefficient of the A term is time-dependent and 
that of the A T term is time- independent. 

4. Both coefficients of the A and A T terms are time- 
dependent. 

The synchronization manifold of the error equation ([2]) 
is locally attracting if the origin of this equation is stable. 
Following the Krasovskii-Lyapunov theory Q , we define 
a continuous, positive-definite Lyapunov functional of the 
form 

1 r° 

V(t) = -A 2 +[i A 2 (t + 9)d9, 1/(0) = (3) 

^ J —T 

where \x is an arbitrary positive parameter, /i > 0. The 
derivative of the functional V(t) along the trajectory of 
the error equation @, 

^ = -(a + K(t))A 2 + bf'(y(t - r))AA T + fiA 2 - //A*, 

(4) 

has to be negative to ensure the stability of the solution 
A = 0. The requirement that ^ < for all A and A r , 
results in the condition for stability as 

a + K(t)> ^f(y(t-T)) 2 +f i = <P(ri. (5) 

Now, 3>(/i) as a function of /i for a given f'(x) has an 
absolute minimum at 

H = {\bf{y{t-T))\)/2, (6) 

with $ min = \bf{y{t-r))\. Since $ > $ mm = \bf'(y(t- 
t))|, from the inequality (|5|), it turns out that a sufficient 
condition for asymptotic stability is 

a + K(t)>\bf'(y(t-T))\. (7) 



It is to be noted that since \i is an arbitrary positive 
parameter due to the definition of the positive definite 
Lyapunov function ([3|), the above stability condition is 
valid only when /i = (\bf'(y(t — r))|)/2 is a constant, 
i. e., only when f'(x) is a constant (in other words 
when the coefficient of A T term in the error equation 
([2]) is time- independent, which corresponds to the cases 
(1) and (3) discussed above). On the other hand if f'(x) 
is time-dependent, then u can be obtained alternatively 
by rewriting Eq. (|5|) as 

b 2 f'(y(t~r)) 2 < ty(a + K(t)-n), (8) 
= -4 [a*- {a + K(t))/2} 2 + {a + K(t)) 2 , 

EE Vfr). 

Now, ^(/i) as a function of [i for a given f'(x) has an 
absolute maximum at 

f i=(a + K(t))/2, (9) 

with ^ m ax = {a + K(t)) 2 . Using this maximum value in 
the right hand side of (8), we obtain the same stability 
condition as that of (O, provided (a + K(t))/2 > since 
[i > 0. Since a > 0, this implies K (t) > —a, that is 
coupling function K(t) should be either positive definite 
or \K (t)\ > a if it is negative. In particular for the case 2, 
since the coefficient of the A term in the error equation 
is time independent (which corresponds to the cases (1) 
and (2) mentioned above), K{t) = k > —a for all t (k : 
const.). 

However, there arises an even more general situation 
where the coefficients of both the A and A r terms are 
time dependent (case 4) , in which case the arbitrary pos- 
itive parameter /i in the Lyapunov functional has to be 
chosen as a positive definite function, u — g(t) > for 
all t. In this case, one has to consider the derivative of 
u = g(t) also in the derivative of V{t) as follows, 

^ = -(a + K(t))A 2 +bf'(y(t-T))AA T (10) 
+ g(t)(A 2 - A 2 T ) + g(t) J A 2 (t + 9)d6<0. 

It is known from the Lyapunov functional that the term 
J_ A 2 (t + 9)d0 is positive and let us suppose that g(t) < 
for all t, then for V(t) < a sufficient condition is that 

- (a + K(t)) A 2 + bf(y(t-r)) AA T (11) 
+g(t)(A 2 -A 2 T )<0, 

- [(a + K{t)) b 2 f'(y(t - r)) 2 /Ag{t) - g(t)] A 2 (12) 
-g(t) [A T - bf(y(t - r))A/2g(t)] 2 < 0. 

The second term in the above equation is negative by 
assumption of g(t) and hence it follows that 

b 2 f(y(t-r)) 2 < Ag(t)(a + K(t)-g(t)), (13) 
= -4lg(t)-(a + K(t))/2} 2 + (a + K(t)) 2 , 

= r( 5 (t)). 
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FIG. 1: (Color online) The time trajectory plot of the vari- 
ables x(t) and y(t) of the coupled time-delay systems {TJ in- 
dicating complete synchronization between them, (a) Piece- 
wise linear time-delay system, ([TJ and (|14p . for the parame- 
ters a — 1.0, b — 1.2, r = 25.0 and for the constant coupling 
fei = &2 = 0.9. (b) Ikeda time-delay system, JTJ and (|15[1 . 
for the parameters a — 1.0,6 — 5,t — 2.0 along with the 
constant coupling k\ = = 5.0. (c) Piecewise linear time- 
delay system for the same values of the system parameters as 
in Fig. [1J, with the square wave coupling rate ki = 0.9 and 
&2 = 1-0. (d) Ikeda system for the same values of the system 
parameters and with the square wave coupling rate ki = 5.0 
and &2 = 6.0. 



Consequently we obtain the same stability condition as 
in Eq. ([7]) with the maximum of T, T max = (a + Kit)) 2 , 
occuring at g(t) — (a + K(t))/2 > 0, along with the 
condition g(t) < 0, that is K(t) = dK(t)/dt < 0, for all 
t. 

Note that our above analysis holds good in case (4), 
only for K(t) < and is not valid for K(t) > 0. For the 
latter case, we are unable to obtain a sufficiency condition 
yet. Consequently the cases (l)-(3) cannot be treated as 
special cases of the most general case (4) at present. 

Thus, we have shown that the same general stability 
condition, Eq. ([7]) , is valid for all the four cases that arise 
in the error equation @ corresponding to the synchro- 
nization manifold of the unidirectionally coupled time- 
delay systems with a restriction in case (4). 

In this section, we will provide numerical confirmation 
of the above stability analysis for all the four cases using 
appropriate nonlinear functional forms f{x) and suitable 
coupling K(t) in the coupled time-delay systems ([T]). For 
this purpose we will consider the nonlinear functions f{x) 
as the piecewise linear function, which has been studied 
in detail recently 0, EE QJ] , 



o, 




x < -4/3 


1.5x - 


2, 


-4/3 < x < -0 


X, 




-0.8 < x < O.t 


1.5x + 


2, 


0.8 < x < 4/3 


o, 




x > 4/3, 



and 

f(x) = sin(x(t-T)), (15) 



which is the paradigmatic Ikeda model [6J] . We have fixed 
the parameters as a — 1.0, b — 1.2 and r = 25.0 for the 
coupled piecewise linear time-delay system defined by ^ 
and (H3J), for which the uncoupled systems exhibit a hy- 
perchaotic behavior with nine positive Lyapunov expo- 
nents [IEEE EE El- For the coupled Ikeda systems ^ 
and (JT5J) , the parameters are chosen as a = 1.0,6 = 5.0 
and t = 2.0 where the uncoupled individual Ikeda sys- 
tems exhibit a hyperchaotic behavior with three positive 
Lyapunov exponents (T^. 

We choose the coupling function K(t) as a square wave 
function represented as [201 ] 

K(t) = {{t Q , ki), (ti, hi), (t 2 , h), (t a , fc 2 ), . . .}, (16) 

where tj = to + (j — l)r s ,j > 1 is the switching instant, 
k\ > 0, ki > with k\ ^ ki. For constant coupling, 
K(t) = k\ = ki- On the other hand, either h\ = 
or ki = 0, then the coupling is called an intermittent 
coupling/control which is now being widely studied in 
the literature [2lj . 

First, we use the piecewise linear function (fl"4"|) , and 
the constant coupling K(t) = hi = fe- It is clear from 
the form of the nonlinear function f(x) and the cou- 
pling that both coefficients of the A and A T terms in 
the error equation, ©, are constant (case 1) and conse- 
quently fi can either be chosen as fj, — (\bf'(y(t — t))\)/2 
or fi = (\a+K(t)\)/2. The time trajectory of the variables 
x(t) and y(t) of the coupled piecewise linear time-delay 
systems, ([1]) and (fl4"|) . are shown in Fig. QJi indicating 
complete synchronization between them for the coupling 
strength k = k\ = ki = 0.9 satisfying the stability condi- 
tion a + k > bf'(y(t — t)) = 1.56. Here, the other system 
parameters are fixed as noted above. 

Next, we analyse the function f(x) = s'm(x(t — r)), 
given by (|15J) . of the Ikeda system with constant coupling, 
which corresponds to the case 2 where the coefficient of 
the A T term in the error equation is time-dependent, 
while that of the A term is still time-independent and 
hence [i can take the form /i = (a + K(t))/2 with K(t) > 
0. The coupling strength is fixed as k = k\ — ki = 5.0 
such that the stability condition a + k > bf'(y(t — r)) = 6 
is satisfied. The variables x(t) and y(t) of the coupled 
Ikeda systems, ([1]) and (IT5|) . are plotted as a function of 
time in Fig. [TJd demonstrating complete synchronization 
between them. 

Again, we consider the piecewise linear function (|14[) . 
and the same parameter values as in the case 1 but with 
the square wave coupling K(t) chosen as k\ = 0.9 and 
ki = 1.0 such that the stability condition is satisfied 
for all t. The switching instant t s between k\ and ki for 
the square wave coupling rate is fixed as r s = 1.0. This 
situation corresponds to the case 3, where the coefficient 
of the A term in the error equation is time-dependent, 
while that of the A T term is time-independent and as 
a result \i can be fixed as fi = (\bf'(y(t — r))|)/2. The 
time trajectory of the variables x(t) and y(t) are shown in 
Fig. [TJj indicating complete synchronization. Note that 
here K(t) > and the stability condition (UJ is indeed 
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satisfied. 

Finally for the more general case where both coef- 
ficients of the A and A r term of the error equation 
are time-dependent, fi = g(t) can be given as g(t) = 
(a + K(t))/2 for the chosen form of the square wave cou- 
pling K(t) with K(t) > and K(t) = 0. Figure EH is 
plotted for the same values of the system parameters as 
in Fig. [TJd with k\ = 5.0, &2 = 6.0 and t s — 1.0 satisfying 
the stability condition ([7]), indicating complete synchro- 
nization between the variables x{t) and y{t). 

Asymptotic stability of synchronized state in a uni- 
directionally coupled time-delay system is studied using 
the Krasovskii-Lyapunov theory. We have shown that 
the same stability condition is valid for all the four cases 
that arises due to the nature of the coefficients of the A 



and A T terms in the error equation corresponding to the 
synchronization manifold. In particular, we have shown 
that the same general stability condition is valid even for 
the general case where both coefficients of the A and A T 
terms in the error equation are time-dependent, which 
is of high importance for various applications. We have 
also numerically confirmed these results using appropri- 
ate examples along with suitable coupling configuration. 
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